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Periodic and Quasi-Periodic Compensation 
Strategies of Extreme Outages caused by 
Polarization Mode Dispersion and Amplifier Noise 

Vladimir Chernyak, Michael Chertkov, Igor Kolokolov, and Vladimir Lebedev 



Abstract — Effect of birefringent disorder on the Bit Error 
Rate (BER) in an optical fiber telecommunication system subject 
to amplifier noise may lead to extreme outages, related to 
anomalously large values of BER. We analyze the Probability 
Distribution Function (PDF) of BER for various strategies of 
Polarization Mode Dispersion (PMD) compensation. A compensa- 
tion method is proposed that is capable of more efficient extreme 
outages suppression, which leads to substantial improvement of 
the fiber system performance. 

Index Terms — Optical Fiber Telecommunication systems, 
PMD-compensation, Polarization Mode Dispersion (PMD), Bit- 
Error-Rate (BER), Probability Distribution Function (PDF) 

Polarization Mode Dispersion (PMD) is an essential impair- 
ment for modern optical fiber systems [1], [2], [3]. There- 
fore, dynamical PMD compensation became an important 
subject in modern communication technology [4], [5], [6], 
[7]. Optical noise generated in optical amplifiers represents 
another impairment that may not be reduced/compensated and, 
therefore, should be also considered in any evaluation of a 
fiber system performance [8]. BER calculated for a given 
realization of birefringent disorder by means of averaging 
over the amplifier noise statistics constitutes an appropriate 
object to characterize joint effect of the two impairments. In 
two preceding papers [9], [10] we have demonstrated that the 
probability of extreme outages (values of BER much higher 
than typical) is substantially larger than one could expect 
from naive Gaussian estimates singling out effects of either 
of the two impairments. The natural object of interest is the 
PDF of BER and, specifically, the PDF tail corresponding 
to anomalously large BER. In [9] we have developed a 
consistent theoretical approach to calculating this tail. The 
case when no compensation is applied and also the effect of 
the simplest "setting the clock" compensation on the PDF tail 
suppression have been discussed in [9]. Then our investigation 
was extended to study effects of the standard first- and higher- 
order compensations on extreme outages [10]. In the present 
letter we propose a compensation scheme that appears to be 
more efficient in reducing the extreme outages compared to 
the traditional high-order compensation scheme with the same 
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number of the compensating degrees of freedom. 

We consider the return-to-zero (RZ) modulation format, 
when optical pulses are well separated in time t, and thus 
can be analyzed as individual objects. We represent the pulse 
intensity measured at the system output as 



AtG{t) \K&(t)\' 



(1) 



where G{t) is a convolution of the electrical (current) fil- 
ter function with the sampling window function. The two- 
component complex field \P(t) describes the output optical 
signal (the components correspond to two polarization states of 
the signal). The linear operator JC in Eq. {3 represents optical 
filtering, it may also account for a compensating device. The 
compensating part of the linear operator, JC C , is applied first, 
i.e. before filtering described by JCf, resulting in JC = JCfJC c . 
Ideally, / takes two distinct values depending on whether the 
information slot is vacant or filled. However, the impairments 
enforce deviations of / from those fixed values. If the output 
signal intensity exceeds the decision level Id, then "1" is 
associated with the slot, otherwise the slot is labeled by 
"0". Sometimes the information is lost, i.e. the initial "1" is 
detected as "0" at the output or vise versa. The BER is the 
probability of such events that naturally depends on a specific 
realization of birefringent disorder in the fiber. BER must be 
extremely small to guarantee successful system performance. 
It has been demonstrated in [9] that anomalously high BER 
originates solely from the "1 — * 0" events. We denote the 
probability of such events by B and study its sensitivity with 
respect to disorder. 

In this letter we restrict ourselves to the linear propagation 
regime, when the output signal \P(t) can be represented 
as a sum of two contributions: ip, related to the noiseless 
evolution of the initial pulse, and the noise-induced part <fi. 
We consider the cases of distributed or, alternatively, lumped 
amplification with the fiber length Z essentially exceeding the 
inter-amplifier separation (span length) within the same frame- 
work. <p becomes a zero-mean Gaussian variable, completely 
characterized by its two-point correlation function 



(0a(<l)^(<2)) = DeZSapSfa ~ is), 



(2) 



insensitive to a particular realizations of birefriengent disorder 
and chromatic dispersion in the fiber. The product D^Z 
represents the amplified spontaneous emission (ASE) spectral 
density accumulated along the fiber. The coefficient is 
introduced into Eq. (0 to reveal the linear growth of the ASE 
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factor with Z [8]. The noise-independent part of the signal is 
governed by 



d z ip — m(z)d t ip — id(z)d 2 <p = 0, 



(3) 



z and d being the coordinate along the fiber and chromatic 
dispersion. The birefringence matrix can be represented as 
rh — hj&j, where hj is a real three-component field and 
<jj are the Pauli matrices. Averaging over many states of the 
birefriengent disorder, any fiber is going through over time, 
or (alternatively) over the states of birefringence of different 
fibers, one finds that hj(z) is a zero-mean Gaussian field 
described by 



(hi(zi)hj(z 2 )) = D rn dij8(zi - z 2 ). 



(4) 



If birefringent disorder is weak the integral H = j Q dz h{z) 
coincides with the PMD vector. Thus, D rn — fc 2 /12, where k 
is the so-called PMD coefficient. 

In an operable communication system, typical damage 
caused by disorder and noise must be small, i.e. typically both 
impairments can cause only a small distortion to a pulse, thus, 
the optical signal-to-noise ratio (OSNR) and the ratio of the 
squared pulse width to the mean squared value of the PMD 
vector are both large. OSRN can be estimated as Iq/(D^Z) 
where Io = / di |<Jo(^)| 2 is the initial pulse intensity, the 
integration being performed over a single slot populated by an 
ideal pulse, encoding "1". Typically, B fluctuates around B , 
the zero-disorder (hj — 0) value of B. A convenient auxiliary 
dimensionless object, T = (D^Z) \n(B/B n )/I , depends on 
the birefringent disorder, the the initial signal shape, as well as 
the details of the compensation and detection procedures, it is, 
however, insensitive to the noise. Since OSRN is large, even 
weak disorder can generate strong increase in the value of B. 
This is why a perturbative (with respect to h) calculation of 
r gives the most essential part of the PDF S(B) of B. If no 
compensation is applied, one gets V ~ H^/b, b being the pulse 
width, and the initial signal is assumed to be linearly polarized. 
In the simplest case of the "setting the clock" compensation 
one arrives at T ~ (H 2 + H%)/b 2 . This yields the power-like 
tail of the PDF of B [9]. Higher-order compensation leads 
to T ~ (H/b) p , where p is an integer exceeding by one 
the compensation degree provided no additional cancellations 
occur, and one gets the following asymptotic expression (tail) 
for the PDF S(B) of B [10]: 

kiS = -ti p b 2 [D^Zln(B/B Q )/I } 2 / p /[D m Z], (5) 

where /i p is a dimensional coefficient. Therefore, as an- 
ticipated, compensation suppresses the PDF tail. However, 
applying high-order compensation is not too efficient, since 
the decrease of S is mild as p grows. 

The main purpose of this letter is introducing more efficient 
compensation strategies with the same number of compen- 
sating degrees of freedom. As a first example consider the 
following "periodic" scheme. One divides the optical line 
into N segments, each of the length I = Z/N, and apply 
the first-order compensation at the end of each segment (as 
schematically shown in the upper panel of Fig. [0 with "c" 
denoting the compensating elements). The noise-independent 
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Fig. 1. Cartoon scheme of fiber- line elements installation correspondent to 
the periodic and quasi-periodic compensation strategies. 



part of the compensated signal for the "periodic compensation" 
strategy is determined by 



JC c <f = exp(ir]dt)lC 1N U N . . . KuU-^^t) 
fill 



U n = Texp 
/Ci„ = exp 



(n-X)l 

i 



dz hj(z)cjjdt 



(n-l)l 



dz hj(z)(jjdt 



(6) 
(7) 

(8) 



where 3*b(i) is the input signal profile, r\ = J Q dzd(z) is the 
integral chromatic dispersion, and the ordered product on the 
rhs of Eq. is taken over all the N segments (T cxp is the 
standard notation for the so-called ordered exponential). The 
exponential factor /Ci„ represents the first-order compensation 
at the end of the n-th segment. 

This "periodic" compensation is not too convenient since 
it requires installation of compensating elements at multiple 
places along the fiber. However, one can naturally modify this 
scheme and have the same compensating elements inserted 
subsequently but all at once at the fiber output as it is shown 
in the lower panel of Fig.^ If tne disorder profile hj is known 
(technically such end-point measurements are possible through 
the anti-Stokes refraction technique [11]) one can have an end- 
point, but multiple, compensation as JC C = n^-ii' leading to 
the following "quasi-periodic" modification of Eq. 0: 



!C c <p = exp(i?79 t 2 )/Cii . ..JC 1N U N . . . Ux& (t) . 



(9) 



The natural idea behind this "quasi-periodic" compensation 
is obvious: to construct (in the compensating part) the best 
possible approximation (with the given number of the com- 
pensating degrees of freedom) for the inverse of the ordered 
exponential Un ■ ■ ■ U\- 

Note that the (quasi) periodic compensation does not in- 
fluence the noise-dependent part of the signal, i.e. K, c cf> has 
the same correlation function (|2) as <p. Therefore, one arrives 
at the same expression \u(B/Bq) = TIq/ (D^Z), with a new 
fa-dependent factor T. Furthermore, in the region of our main 
interest T can be analyzed perturbatively, just as in [9], [10]. 
Expanding the factors in Eq. (|6} up to the second order and 
making use of Eqs. (II 121 one derives: 



b 2 ^ 



ii 
dz 



dz' [h^z^iz') - h 2 (z)h 1 (z% (10) 
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where a n = (n — 1)1. Here, the dimensionless coefficient \J 2 
is related to the output signal chirp produced by an initial 
chirp and the nonzero integral chromatic dispersion rj. As 
follows from Eq. l|9}, the same expression dlOi is obtained 
in the second order for the quasi-periodic case. Substituting 
Eq. JlOi into the expression for B and evaluating the PDF 
of B, with the Gaussian statistics of h described by Eq. 10), 
leads to the following expression for the tail of the PDF of B: 



S(B) dB 



BgdB 



(11) 



Eq. (O holds for \n(B/B ) > ii' 2 D m I /[D(:b 2 ]. The expo- 
nent a in Eq. il Q contains an additional factor N compared 
to the expression for the first order end-point compensation, 
i.e. the (quasi) periodic compensation makes the tail of S(B) 
steeper. It is instructive to compare the outage probability 
for the periodic case with the case of higher-order end-point 
compensation described by Eq. One finds that for higher- 
order compensation, i.e. when N ~ p! 2 ln(S*/ B )I / (D^Z), 
the (quasi) periodic scheme becomes more efficient compared 
to the straight 7V-th order compensation scheme. If the output 
signal is not chirped, [i! 2 = and the leading term in the 
expansion of V in h/b is of the third order. Additional filtering 
efforts can be made to enforce the output pulse symmetry un- 
der the t — ► —t transformation, thus removing the third-order 
term. Then the leading term in T will be of the fourth order in 
h/b. Finally, even better compensation can be achieved if the 
standard high-order compensation approach and the (quasi) 
periodic ones are combined, i.e. if in the (quasi) periodic 
setting, considered above, one uses higher-order compensation 
instead of the first order one. Formally, this hybrid case means 
that the first-order compensation operators JC\ n in Eqs. ( I6I9> 
should be substituted by higher-order compensation operators 
K, cn . In the hybrid periodic case T can be written as a sum 
of r„ ~ (JJ 1 dz h/b) p , and, since h is short-correlated, T n 
related to different segments are statistically independent. This 
leads to the following expression for the PDF tail 



In S(B) 
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valid at D^Z J I \n{B / B ) > N 1 -P/ 2 (D m Z/[fi p b 2 ])P/ 2 . 

Note, that an important computational step, leading to our 
major results in Eqs. (II 1I12> . was evaluation of T pertubatively 
in h. Besides, in the periodic case T is a direct sum of 
each segment contribution T n , and the perturbative treatment 
applies separately to each T„, requiring the weakness of 
the PMD effect at each segment only, i.e. D rn Z/N < b 2 . 
Therefore, one concludes that even an optical line with not 
really operable (without compensation) characteristics (D m Z 
which is of the order or larger than b 2 ) can still be used for 
transmission if N is sufficiently large. Moreover, this obser- 
vation on the applicability of Eqs. il 111 2t also extends to the 
quasi-periodic case, in the sense that Eqs. (Ill I12i provide an 
upper bound for the PDF of BER. This is due to an additional, 
oscillatory with h, suppression of T n in the quasi-periodic 
case vs periodic. This suppression is especially important for 
segments strongly separated from their compensating counter- 
segments. 



For illustration purposes let us briefly discuss an example 
of a fiber line with typical bit error probability, _B = 10~ 12 , 
and p! 2 — 0.14. Assume also that the PMD coefficient, 
k = ^/12D m , is 1.5 ps/Vhm, the pulse width, 6 = 25 ps, 
and fiber length Z — 2, 500 km. Then the dimensionless 
parameter, D m Z/b 2 measuring the relevant strength of the 
PMD effect is 0(1), i.e. without compensation PMD effect 
is large, pulses are destroyed and no successful transmission 
is possible: S(B) = 0(1) for any B > Bo. If, however, the 
(quasi) periodic compensation with N = 10 compensation 
units is utilitized the relevant strength of the PMD effect is 
essentially reduced, so that D m Z/[b 2 N] ps 0.1, and S(B) 
starts to decay with B at B > Bq. The system performance 
can be evaluated in terms of the outage probability C,the 
probability for B to be larger than O = Jg dBS(B). 
One derives from Eq. (lilt that in the (quasi) periodic case, 
O « 0.06 for B* = 10~ 8 and O « 0.02 for B* = 10~ 6 , i.e. 
the system performance is essentially improved (to become 
not yet perfect but already satisfactory). 

To conclude, in this letter we have proposed a (quasi) 
periodic compensation scheme which appears to be a strong 
alternative to the standard higher-order compensations. The 
efficiency of the scheme has been demonstrated. Even though 
technical implementation of this procedure needs an expensive 
equipment, we anticipate that if this compensation technique 
is implemented the reduction in the probability of extreme 
outages will guarantee an essential overall benefit. 

The authors thank I. Gabitov for numerous valuable discus- 
sions. 
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